Out of equilibrium transport through an Anderson impurity: Probing scaling laws within the 

equation of motion approach 
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We study non-equilibrium electron transport through a quantum impurity coupled to metallic leads using the 
equation of motion technique at finite temperature T. Assuming that the interactions are taking place solely 
in the impurity and focusing in the infinite Hubbard limit, we compute the out of equilibrium density of states 
and the differential conductance G2 (T, V) to test several scaling laws. We find that G2 {T, V) /G2 {T, 0) is 
a universal function of both cV/Tk and T/Tk, being Tk the Kondo temperature. The effect of an in plane 
magnetic field on the splitting of the zero bias anomaly in the differential conductance is also analyzed. For a 
Zeeman splitting A, the computed differential conductance peak splitting depends only on A/Tk, and for large 
fields approaches the value of 2A. Besides the traditional two leads setup, we also consider other configurations 
that mimics recent experiments, namely, an impurity embedded in a mesoscopic wire and the presence of a third 
weakly coupled lead. In these cases, a double peak structure of the Kondo resonance is clearly obtained in the 
differential conductance while the amplitude of the highest peak is shown to decrease as \n{eV/TK)- Several 
features of these results are in qualitative agreement with recent experimental observations reported on quantum 
dots. 

PACS numbers: 72.15.Qm,73.21.-b,73.23.Hk 



I. INTRODUCTION 

The Kondo effect has been experimentally tested in sev- 
eral devices, ranging from atoms and molecules in contact 
with metallic reservoirs to quantum dots embedded in com- 
plex mesoscopic structures.'"^ This many body effect occurs 
when the ground state of the atom, molecule or quantum dot 
(from hereon 'the impurity') is degenerated — typically due to 
spin degeneracy. While at high temperature the impurity be- 
haves as a free spin, coiTelations build up as the temperature T 
is lowered. As a consequence, below a characteristic temper- 
ature Tk, called the Kondo temperature, the impurity's spin is 
screened by the spins of the electrons in the host structure. As 
a result, the Fermi liquid picture is recovered as T ^ 0. For 
a system in equilibrium, a characteristic feature of the Kondo 
regime is the emergence of a sharp resonance in the impurity 
spectral density at the Fermi level. In the absence of exter- 
nal magnetic fields, Tk is the only low energy scale and all 
the (zero bias) transport properties obey a universal scaling as 
a function of T/Tk- This effect, which is theoretically well 
understood, has been reported in a variety of experiments. "*"™ 

External magnetic fields B or bias voltages V introduce 
new energy scales that could be of the order of Tk- In 
these conditions, several recent works focused on the search 
of scaling laws to test whether the non-equilibrium Kondo 
physics remains universal. In particular, measurements of 
the two terminal differential conductance G2 {T, V) in sin- 
gle channel quantum dots seem to be well described, at rel- 
atively low V and low T, by a universal scaling law char- 
acterized by two parameters.^ The universal character of the 
magnetic field induced splitting of the zero bias anomaly in 
the differential conductance, has also been tested in a series 
of experiments.'"'"''' At low B, the splitting decreases with 
Tk and the data are compatible with the predicted univer- 
sal behavior, being B/Tk the scaling variablefi^iS A recent 



publicationr^ however, reported that at high B there is a 
crossover regime in which the behavior of the splitting on Tk 
seems to be inconsistent with the universal character 

From the theoretical point of view, the main ingredients re- 
quired to compute the electronic cuiTent and related quantities 
are the retarded and lesser Green functions^" of an interact- 
ing quantum impurity embedded in a metallic networkii2i22 
Much effort has been devoted in recent years to find reli- 
able approximations for these functions employing different 
techniques. Most of the approaches consider the single im- 
purity Anderson model, in which a single correlated elec- 
tronic state is connected to two leads. For this configura- 
tion, some theoretical considerations predict an universal scal- 
ing of G2{T, V) in temperature and applied bias but disagree 
about the coefficients and parameters that determine the scal- 
ing function. '^•^^•^"' Other studies that employ non-equilibrium 
renormalized perturbation theory, need to include valence 
fluctuations in order to reproduce the simple scaling formula 
for G2{T, V).^^ The non-crossing approximation (NCA)2Ii^ 
predicts the splitting of the out of equilibrium Kondo reso- 
nance into two peaks, but fails to reproduce other aspects of 
the Kondo correlations and the spectral density presents an un- 
physical peak at the Fermi level in the presence of a magnetic 
field. Discrepancies in the results may be due to the fact that 
these approaches, although well tested at equilibrium, have 
several weak points when extended to the out of equilibrium 



situation 
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Using the equation of motion (EOM) technique with the de- 
coupling scheme proposed by Lacroix,^*^ the authors of Refs. 
131.32.1 have recently obtained an explicit analytical expres- 
sion for the impurity's retarded Green function in a general 
case in which the impurity is embedded in a mesoscopic struc- 
ture at equilibrium. These studies showed that it is crucial to 
include in the EOM all the correlations emerging from the 
truncated scheme in order to reproduce the low temperature 
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FIG. 1: (Color online) a) Schematics of the two terminal quantum 
dot setup; b) quantum dot coupled to a biased ballistic wire and a 
non-invasive probe; c) quantum dot in a three terminal setup. 



Kondo correlations that other approximations like the NCA 
fail to reproduce. 

The purpose of the present work is to extend this approach 
to the out of equilibrium situation to analyze the transport 
properties of a single quantum impurity and test some rele- 
vant scaling laws. First, we shall revisit the conventional two 
leads setup displayed in FiglT^) and show that the results for 
G2{T, V) are indeed universal; they depend only on cV/Tk 
and T/Tk- We also scale our results with the simple empiri- 
cal (and approximated) law recently used in the literature^ to 
scale the experimental data corresponding to a single channel 
quantum dot and find a similar behavior. In addition, we dis- 
cuss the effect of an in-plane magnetic field B on the splitting 
of the zero bias anomaly in the differential conductance. In 
agreement with Ref. 116(1 , and in contrast to Ref. |19], we 
find that the magnitude of the splitting, 25V, is limited by the 
Zeeman energy, 5V < gjisB. 

For single level impurities G2 (T, V) is essentially deter- 
mined by the impurity's spectral density Pd{<^)- However, 
the conductance does not provide enough information to ex- 
tract the V dependence of Pd{'-^)- For this purpose, alterna- 
tive set ups that include an additional lead have been theo- 
retically proposed. ^^^^ Conductance through the additional 
weakly coupled lead provides a direct probe of the out of equi- 
librium spectral density in the Kondo regime. Closely related 
experimental setups have been designed in order to achieve 
that. In Ref. iQ] a quantum dot was coupled to a mesoscopic 
wire formed in a 2DEG that was brought out of equilibrium by 
an applied bias voltage. More recently, in Ref. |4] the Kondo 
effect was measured in a ring-shaped dot that was connected 
to two strongly and one weakly coupled leads. Both experi- 
ments demonstrated unambiguously the splitting of the zero 
bias conductance peak, attributed to the splitting of the Kondo 
spectral density due to the double step distribution function of 
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FIG. 2: (Color online) Equilibrium spectral density for Ed = — 4r, 
r = 0.015i3 y r = Tif. Inset: Detailed behavior of the Kondo peak 
for different temperature values, T = 0.1, 1, IQTk 

the caiTiers in the leadsi^ 

We shall discuss two related setups (Figs, ^p) and[T];)) that 
allow the detection of the out of equilibrium spectral density 
of the quantum impurity: i) in FigH})) we mimics the set up 
used in Ref. [2] where a quantum dot is coupled to an out of 
equilibrium (biased) ballistic wire and a third contact is used 
as a weakly coupled (non-invasive) probe; ii) in Fig[T]b) an 
extra lead has been added to the two lead configuration repro- 
ducing the setup of Ref. [4] . 

The paper is organized as follows: Section HIl starts with a 
description of the model and the EOM approximation. In sec- 
tions Hn] and HV] we present the results for the two and three 
terminal configurations, respectively. Section [V] contains a 
summary and conclusions. An appendix with some details 
of the out of equilibrium calculation is also included. 

11. THE MODEL AND GENERAL FORMULATION 

The Hamiltonian for a single level impurity coupled to a 
mesoscopic structure as the ones depicted in Fig[T] is given 
by iJ = Hd + Hr- The first term in H describes the isolated 
impurity, 

Hd='^ £dcrnda + U UdlUdi , (1) 

Udcr = d\d^ is the occupation number operator and dj. cre- 
ates an electron with spin projection a at the impurity level 
with energy £da — £d ^ 0'9PbB/2 where Ed is the single 
particle energy in the absence of the external magnetic field 
B. The last term in Eq. ([T]) describes the Coulomb interac- 
tion U. In the following we take the limit U ^ 00, which 
forbids the doubly occupancy of the dot level. The other term 
Hr describes the non-interacting leads and its coupling to the 
impurity, 

Hr = ekacika^^aka + (14fccJ,fc„d<x + h-C), (2) 

where cj^^,^ creates an electron with momentum k and spin a 
in the a-lead. The leads are in contact to reservoirs at differ- 
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ent chemical potentials that bring the system out of equi- 
librium. In the case of the two lead configuration depicted in 
Figdia), the label a stands for the left (L) and right (R) con- 
tacts with chemical potentials ji^ and fin, respectively. In this 
configuration, the current / flowing through the impurity can 
be expressed in terms of the impurity's spectral density^^ 



= - - VlmG5^(w) 

IT 



(3) 



where G'[i^{io) is the retarded Green function of the impurity. 
Using the EOM technique (see the Appendix for the details) 
we obtain the following expression for the non-equilibrium 
(retarded/advanced) Green function, 



l-{n,)-P;'\0Js) 



(4) 



with Ljg = u + agfisB (a — — cr). 



(5) 



and 



a/r 

w ± lO 



■Auj'. 



(6) 



Here, — T^^k l^afcP ^ ^i'^ ~ ^ka) defines the cou- 
pling to the leads for which we assume a constant density 
of states of 1/2D per spin, F = J^a^o' '^^'^ fa{'^) = 
(exp — /iQ,)] + 1)^^ is the Fermi distribution function of 
the a-lead with (3 — l/ksT. Notice that the set of equations 
and (|6]l, which have to be solved self-consistently using 
numerical methods, have the same form as in the equilibrium 
cas©^ but with an effective distribution function 



(7) 



The reason for this is that, within the decoupling scheme 
used for solving the EOM, the impurity's lesser Green func- 
tion can only take the pseudo-equilibrium form G'^^{uj) = 
/(w) [GaA^) - G5^„(w)]— see the appendix. 

At equilibrium (/(w) = /(w)), the EOM approach leads 
to the impurity density of states Pcq(i^) shown in Fig. |2] It 
presents a clear Kondo peak near (ui = 0). A detail of the 
Kondo peak for different temperatures is shown in the in- 
set. As stated in Refs. ll3lll32ll . as T the impurity's 
spectral density satisfies the Fermi liquid unitary condition, 
Pcq{io = 0) = 2/7rr. Here, however, that limit is reached 
only logarithmically as the exact low temperature depen- 
dence is not recovered within this approach.^' 

In what follows we use these results to probe the scal- 
ing laws of the transport properties in different configurations 
with and without an external magnetic field B. 
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FIG. 3: (Color online) Linear conductance as a function of T/Tk 
for different values of Sd and for F = 0.015-D. There is a very good 
scaling with the effective Kondo temperature Tk = D exp{-Ked/^)- 
Notice that the scaling is improved as |£d|/r increases and that the 
exact quadratic low T behavior of the conductance is not captured in 
this approach. Inset: Unsealed conductance. 



III. TWO TERMINAL CONFIGURATION AND SCALING 
LAWS 



The current / flowing through the impurity in the two ter- 
minal setup can be written as^- 

^ = f f^Tf^ / '^'^(^)[/^('^) - • (8) 

The two terminal linear conductance, given in terms of the 
equilibrium density of states pcq{uj), is then 



G2(T,0) 



/(-f)pe,Hd.. (9) 



The results for G2(r, 0) vs T/Tk for four different values 
of Ed/F and Ti^ = Vr are shown in Fig. [3] Here, Tk = 

exp (TTErf/F) is the relevant energy scale, or the effective 
Kondo temperature, that is obtained within this approach.— 
Unsealed data is shown in the inset for comparison. For all 
the parameters studied, we obtain an excellent scaling in the 
temperature range Q.ITk < T < IOTr-. Note that for large 
|£d|/r the scaling extends over a wider range of temperature. 
In addition, it is apparent from the figure the absence of the 
characteristic Fermi liquid behavior expected at low T. 
Instead, a logarithmic behavior is obtained. 

The non-equilibrium density of states pd{'^) with /i^ = 
—jiR — eV/2 is shown in Fig. |4]for several values of the 
bias voltage V and T = Tk- As V increases the Kondo peak 
splits with two maxima at approximately Wp ~ ±ey/2. This 
is shown in detail in the inset of the figure, where the dashed 
lines, corresponding to the expression uj-p — ±eV/2 + uj^, 
are a guide to the eyes. As the Kondo peak has an asymmet- 
ric lineshape, due to the background contribution to the spec- 
tral density, there is a difference in the relative height of the 
two non-equilibrium peaks. It is worth to point out that their 
heights will also change in the case of asymmetry coupling. 

The differential conductance G2{T,V) = dl/dV is ob- 
tained as the numerical derivative of Eq. The result, nor- 
malized by G2(T, 0), is plotted in Fig. |5^) as a function of 
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FIG. 4: (Color online) Non-equilibrium spectral density. The figure 
shows the splitting of the Kondo peak for different values of the bias 
voltage, eV G [0, 2QTk\- Here Sd = -4r and V = 0.015D. Inset: 
Kondo peak position ojp as a function of the bias. The dashed lines 
correspond to uj-p = ±eV/2 + uj"^ where lu^ is a fitting parameter. 



(ey/T/f)^ for different temperatures and five different values 
of |£d|/r. All the data corresponding to the same T /Tk col- 
lapse into a single curve. This shows that the EOM approach 
is able to capture the universal scaling as a function of cV/Tk, 
as suggested in Ref Isill . Furthermore, we are able to fit the 
numerical data by the following expression 



G2{T,V) 
G2(T,0) 



eV_ 
TV 



(10) 



where the coefficients a„ are a function of T/Tk only. 
This is shown in Fig. |3b), where we have plotted [1 — 
G2(T, V)/G2{T, as a function of {eV/TkY where i = 
[ai+a2{eV/TKf+az{eV/TKY + ai{eV/TKf]. The good 
scaling shows that, within the EOM approach. 



G2(T,0) 



T eV 



(11) 



is indeed a universal function of both cV/Tk andT/Tx- This 
is valid only in the U oo limit , as it is known that for finite 
U the function has non-universal contributions — it depends 
on the asymmetry of the couplings, for instance. ^''•^'^ Finding 
the analytical form of T is, however, difficult. 

The scaling properties of G2{T,V) were recently exam- 
inated by experimental groups in both quantum dots^ and 
molecular systems.— It was shown that the differential conduc- 
tance approximately follows a rather simple universal scal- 
ing relation, empirically chosen so that the low T and low 
V regime agrees with the quadratic behavior expected from 
Fermi liquid theory. Namely, 



1 - 



G2iT,V) 
G2(T,0) 



ay 



eV_ 
TV 



(12) 



where ay — CTa[l + cxi'^/a — 1){T/Tx)'^] ^ — here ay 
plays the role of ai in Eq. ( fTOt . The parameters ct, 7 and 
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FIG. 5: (Color online) a) Normalized differential conductance as a 
function of (eV/TK^ for T/Tk G [0.25, 2.25] and ea/r = -14/3 
triangles), — 16/3 (circles),— 6, 20/3 and 7.07 (lines). Open symbols 
correspond to T < Tk and the arrow indicates the direction of in- 
creasing temperature. Note that the scaling with bias voltage is rea- 
sonable good even beyond the quadratic regime.; b) Scaling of the 
normalized differential conductance as discussed in the text, c) Scal- 
ing suggested in Ref. |5] , Eq. l iI2t . 



a depend on the nature of the system and experiments have 
found quite different values for quantum dots or molecular 
systemsi^ While in GaAs quantum dots these parameters are 
in reasonable agreement with theoretical calculations, ^"^'S^i^^i^ 
in molecular junctions there are systematic deviations from 
these values whose origin is still unclear."*' Therefore, we use 
the experimental values obtained in Ref.|5fl to scale our nu- 
merical data for the purpose of comparison. The result is 
shown in Fig. EJc). For cV/Tk ^ 0.7 all curves approx- 
imately collapse into a single line. In agreement with the 
experimental observation^ and recent calculations,.^^ devi- 
ations from the [cV/TkY behavior are more pronounced for 
low values of T/Tk- This is related to the fact that the scal- 
ing function ay is too simple to capture the universal scaling 
shown in Fig. |5jb). It is clear then, that a better test of univer- 
sality would be to fit the experimental data with Eq. (fTOl i and 
check the scaling of the coefficients with T/Tk- 

Several experimental groups have recently measured the ef- 
fect of an in-plane magnetic field on the differential conduc- 
tance at finite bias voltages^-i^ii^ with different results. Since 
the Zeeman coupling shifts the impurity levels generating a 
splitting A — gfisB between the two spin states one would 
naively expect a splitting of the Kondo peak in the spectral 
density and a splitting of the zero bias anomaly in the differ- 
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FIG. 6: (Color online) a) Non-equilibrium spectral density for 
A/Tx = 2, eV = A, Sd/V = 16/3 and T/Tk = 1/5 (dashed 
line) and T /Tk ~ 1/10 (solid line); b) Differential conductance 
for different values of the Zeeman splitting, A/Tk £ [0,2] and 
T /Tk = 1/5. The conductance is normalized by its zero bias and 
zero B field value; c) Position of the conductance peak {±5V) as a 
function of the Zeeman splitting for T /Tk = 1/2 (*), 1/4 (A), 1/5 
(•) and 1/10 (■). The open dots are the experimental results of Ref. 
Ilql . Note that the splitting of the conductance peak is well defined 
for A > QATk when T /Tk = 1/10. Inset: Zoom in around the 
experimental data. 



ential conductance. From the theoretical perspective, this is 
indeed the situation for /9oq(a;); it is well established that for 
sufficiently large A the Kondo peak splits in two spin resolved 
peaks at ~ ±A. However, in order to measure the field in- 
duced splitting in a two terminal setup, a finite bias voltage of 
the order of A needs to be applied. As we discussed above, 
this in turns leads to a bias induced splitting and therefore the 
experimental situation is more involved. In fact, in Ref. Ill6ll 
the splitting of the zero bias anomaly was found to be equal to 
2A for large fields while in Refs. IlilUfTl it was found to 
be larger than 2 A. 

Figure |6la) shows the impurity spectral density calculated 
using the EOM approach for a finite bias and a finite Zeeman 
field {eV = A). Two different temperatures and a symmetric 
drop of the bias voltage, fiL — ~IJ-r = eV/2 are considered. 
For the lowest temperature, there are four peaks clearly visible 
at w ~ ±A ± eV/2. The differential conductance 6*2 (T, V), 
normalized to its eV = A = value, is shown in Fig. |6|b) 
for A/Tk e [0,2] and T/Tk = 1/5. For this temperature, 
G2(T, V) presents a local minimum at eT^ = and two max- 
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FIG. 7: (Color online) a) Scanning probe's differential conductance 
as a function of its chemical potential /is for hl ~ eV, hr — with 
eV taking different values in the range [--20Tk , ^OTk] — the curves 
include an offset in the vertical axis for clarity. The position of the 
conductance's peaks and the amplitude of the highest are shown in 
b) and c), respectively, as a function of the bias voltage for Ed = 
-4r(»), -16/3r(B) and -GF (A). The peak amplitude shows a 
\n{eV/TK) dependence. 



ima at ±6V for A/Tk ^ 0.8. The splitting of the zero bias 
anomaly increases monotonically with A being clearly visible 
even though the corresponding peaks in the spectral density 
are barely observable. The monotonic behavior of the con- 
ductance peak position SV as a function of A /Tk is plotted in 
Fig. Wjc) for different temperatures. The experimental data of 
Ref. iflql . which corresponds to T/Tk ~ 1/6, is included for 
comparison. The EOM results corresponding to T/Tk ^ 0.2 
are in qualitative agreement with those data. We note that SV 
is well defined when A is greater than a certain critical value 
Ac— for instance, Ac w OATk for T/Tk = 1/10— which 
increases with T, as expected. We observe that for a large 
Zeeman field 6V A, in agreement with some previous the- 
oretical estimates. We have verified (data not shown) that 6V 
depends only on A/Tk, in contrast to the experimental data 
of Ref. 1 19]. 

It is worth to emphasize that even though G2 {T, V) is not 
simply related to the impurity spectral density, as in the linear 
regime, the value of the observed splitting is already an indi- 
rect evidence of the strong bias dependence of Pd{'-^)- Indeed, 
had we calculated the differential conductance using the equi- 
librium spectral density, Eq. (|9]l, we would have found that 
G2(r ~ 0, 1/) cx /9eq(eF/2) + p^^{-eV/2) and so it would 
present, for large A, two maxima at 5V ^ ±2 A rather than at 
±A. 



IV. THREE TERMINAL CONFIGURATION: PROBING 
THE SPLITTING OF THE KONDO PEAK 

As akeady mentioned, the differential conductance in the 
two lead configuration, G2 (T, V), does not give a detailed in- 
formation of the voltage dependence of the impurity's spectral 
density Pd{^)- To scan the out of equilibrium Kondo reso- 
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nance an extra weakly coupled lead has to be attached to the 
impurity. ^'^"''^ By sweeping the voltage of this extra lead it 
is possible to extract the frequency dependence of the out of 
equilibrium density of states Pd{'^)- Here, we analyze the se- 
tups of Ref. HI illustrated in Figlllb) and (c), respectively. 
The former corresponds to a dot coupled to an out of equi- 
librium wire. An electric current is flowing through the wire 
and a weakly coupled lead connected to the dot (with coupling 
constant Ts) acts as the scanning probe. 

Assuming that the distribution functions of the left and right 
moving carriers in the ballistic wire are given by and 
/i^(cj), respectively, the current flowing from the scanning 
probe through the impurity can be effectively computed em- 
ploying a three lead configuration as the one presented in Ref. 
134]. When the scanning probe is biased to a voltage Vs, the 
current flowing though the impurity is, to lowest order in Fg,^"* 

- ^r, j [f[Lo) - fs{u:)] pdiio) dc-j , (13) 

where /s(w) — l/(exp/3(a; + eVg) + 1) is the Fermi dis- 
tribution of the probe, /(w) is the (out of equilibrium) car- 
rier distribution in the wire and pd is the impurity's spectral 
density in the absence of the probe lead. When the impu- 
rity is coupled with the same coupling constant to left and 
right moving electrons, Vl = F/f , and we have /(cj) = 
{Jl{lu) + fii{u!))/2. The double step distribution /(w) is 
first used to evaluate the spectral density using Eqs.(|4]i and 
(|6]l and then the current Ig. The differential conductance 
G,iT,V)=dIs/dV,is 

GAT,v) ^ -^r.y p.Md^, (14) 

and at low temperatures (T, ~ {2e^ /h)T sPd{eVs). In 
Fig. |2] we present results for the differential conductance 
Gs(T, V) as obtained from Eq. ( fT4l i. We have set = eV 
and PR — 0. The scanning probe conductance as a func- 
tion of the voltage Vs shows a double peak structure due 
to the splitting of the Kondo resonance. As shown in Fig. 
|7lb), the conductance peaks are observed at eVg — pl and 
eVs = PR- The amplitude of the highest peak decreases 
with the wire voltage V and shows a good scaling when plot- 
ted as a function of eV/Tk — in particular for |ed|/F 1. 
Once again, these results are in qualitative agreement with 
the experimental observations.^ However, we obtain a clear 
\n.{\eV\/TK) dependence for the peak's amplitude in the 
range shown in the figure while the experimental data was 
fit to A\n-^{\eV\/TK) + Bln^^ {\eV\/TK).'^ Nevertheless, 
due to the limited range and small values of eV/Tx explore 
in the experiments, the ln(|eV^|/rft-) dependence cannot be 
completely ruled out. 

Similar results are obtained in the case of a single impu- 
rity attached to three leads when one of them is acting as 
a probe, see FigHIc). In that case, f{uj) ~ [TLfhi^) + 
^ rJ r{'^)) I L + Ffl) and thus a left and right asymmetry 
in the couplings changes the relative weights of the two split- 
ted conductance peaks. In particular, for a large asymmetry 
and in the small bias regime, eV < Tk, the unsplitted peak 



tends to be pinned to the chemical potential of the strongest 
coupled lead, as observed in Ref. [4J. 



V. SUMMARY 

We have discussed the out of equilibrium transport prop- 
erties of an Anderson impurity embedded in a multi-terminal 
structure. In order to compute the relevant transport quan- 
tities we have generalized the EOM approach developed in 
Refs. l(3lll32ll to the non-equilibrium situation. Within this 
approach, we have tested the scaling laws recently reported 
in a series of experiments performed in quantum dots in con- 
tact to metalUc leads or wires driven out of equilibrium and/or 
in the presence of an external magnetic field a^i'"^i'^i'^ . We 
found that, within the EOM approach in the U ^ oo limit, 
G2{T,V)/G2{T,0) is a universal function of both bV/Tk 
and T/Tk- Indeed, we have been able to fit the numerical 
data by the expression shown in Eq. ( fTOl i, where the expan- 
sion coefficients a„ are a function only of T/Tk- We then 
conclude that a fit of the experimental data to Eq. ( fTOb instead 
of the simple scaling function ay,^ would be a better test of 
universality. 

In the presence of a Zeeman field we have found a mono- 
tonic behavior of the differential conductance peak position 
5V as a function of A /T k, in agreement with the experimen- 
tal data reported in Ref. lUol . In addition, for large Zeeman 
field we have found that 5V A and that 5V depends only 
on IS./Tk, in contrast to the experimental data of Ref. 1 19]. 

In order to scan the out of equilibrium impurity's spec- 
tral density recently analyzed in a series of experiments^'- we 
have considered, besides the traditional two leads set up, ad- 
ditional configurations, namely an impurity embedded in a 
mesoscopic wire and the three leads configuration in which 
the extra lead acts as a weakly coupled scanning probe. The 
computed scanning probe conductance as a function of the tip 
voltage Vs shows a double peak structure due to the splitting 
of the Kondo resonance. The amplitude of the highest peak 
decreases as a function of the wire's bias with a hi{eV /Tk) 
dependence. In summary, the formalism developed in this 
work provides a reasonable description of the out of equi- 
librium physics of a quantum impurity embedded in a rather 
general mesoscopic setup. More involved setups and possible 
further generalizations are currently under consideration. 
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Note added: After completion of this work we learn of 
two recent manuscripts, Refs. |42.43], that employ a simi- 
lar decoupling scheme of the EOM to compute the impurity 
retarded Green function in the out of equilibrium situation. 
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Appendix: Calculation of GS(tj) using EOM 

In this appendix we present an extension of the EOM results 
of Ref.ll30l-[32ll for the out of equiHbrium case. We use the 
following standard definitions 

{{Ait),B{t'W = -ie{t-t'){A{t)B{t')+B{t')A{t)) 
{{A{t),B{t'))r = i9{t' -t){A{t)B{t')+B{t')A{t)), 

(A.l) 

where the a and r supraindex stand for the advanced and 
retarded Green functions, respectively. For the local Green 



functions in the stationary state we use the notation G^J°'{t — 

t') ^ {{d„{t),dl{t'))Y''' while Gti^) = {{d.,4)y/^ 
indicates its Fourier transform (the same applies to a generic 
green function {{A, B)y/"'). We also define the local lesser 
functions as: G'^^{t — t') — —\{d'l{t')da{t)) and its Fourier 
transform G^^{uj). 



We now consider a local impurity coupled to an arbitrary 
number of reservoirs. Following the standard EOM procedure 
in the infinite U limit, we obtain the following expression for 
the retarded and advanced local Green functions 



+ {{dlrCak,aCpq,a,di)) - {{c\^;gCak,ada , d^))] (A.2) 

I 



where (..) indicates the (out of equilibrium) mean value, 
n„ = c?J.c?cr, a and [3 label the reservoir and gakai'^) = — 
Eafco-)"^ is the reservoir propagator and LUg- — lu + agiisB, 
where B is the external Zeeman field. The decoupling scheme 
proposed by Lacrois^ leads to the following expression 



gj'''{oo)\ G'-J^iu^) = 1 - {n,) - P;'\io,) (A.3) 



with 



(A.4) 

YyJ°'{u}) can be written as the sum of three different contribu- 
tions, = T.\':{u^) + Y^ti^) - Pp''{co)Y.%\u,), 

where 

a.k 

is the non-interacting self-energy and 

Y.l':{^) = Y.V^iJ:^{^)Ue^k.) (A.6) 

cxk 



a.k,f3q 

Here /a(u;) is the Fermi function of lead a, and 

p;/<^(c.) = Y.^o.g'-Jkli^WWk,.) . (A.7) 

ak 



As mentioned above, all mean values need to be evaluated 
in the out of equilibrium situation. In our case, the two rele- 
vant mean values are 



(rfJ-CQfc,o-) = - 

and 



-^y^l'^[9^kaHGM+9lkA^)Gl[M] 

(A.8) 



^a.k,a) ' 



^OikJ3qfct (^aA;o" ) 



'^y-n I ^ g^k^G^Mg^pqA^) + glt,,{u)Gt{^)9l^A^) + glkA^)G^Mg<^A^) 



(A.9) 



where g'^^^i^^) = 2nifa{io)d{uj—eaka) is the non-interacting lesser Green function of the a— lead. Then, the calculation of 
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these mean values requires the knowledge of the impurity's 
lesser green function G^^{uj). An equation for G^^(iuj) could 
be obtained by analytic continuation of the EOM equation if 
the matrix order is preserved. However, in this case it is possi- 
ble to obtain G'^^{lu) by simply imposing that the impurity's 
Green function must satisfy G^^ — C^* where * indicates 
complex conjugation. Indeed, this condition implies that 

(A. 10) 

From the definition of P^{oj), it is possible to show that 



where 



TT Ul — Ul' 



iHaiio) (A. 11) 



(A. 12) 



and 



Here, T„{lo) = tt^^q V^paa{^) = Y.a ^aa{<^) with the 
spin dependent density of states of the a-lead. A similar, but 
cumbersome, expression can be obtained for J7^{uj) — Yj% (w). 
The important point to emphasize is that, as in the case of Eq. 
( lA.l lb . — S° (cl>) is proportional to T-L^- Therefore, 

Hcr(w) = is always a consistent solution of Eq. dA.lOl l. 



Probing that it is the only solution in the general case is not 
straightforward. However, in the absence of an external mag- 
netic field or in the case that the Zeeman splitting is the same 
for the impurity and the band states, Tg{iLia) = ^ai'^), one 
obtains that 



(A. 14) 



Therefore, if we assume that 7^ 0, we arrive to the con- 

clusion that — ImG^g.(aj) = l/rcr(a;), which is clearly an 
unphysical result — note that ?^*(cli) = —"Haii^)- Then, we 
conclude that 7/cr(w) = 0, which corresponds to 



G<,M = -/.M [g:,,h-glm] . 



(A.15) 



This result is nothing but the Ng's ansatz^i which at the this 
level of decoupling of the EOM is the only consistent solution 
for G^^{ll>). Inserting Eq. dA.lSI l into Eqs. dA.St y iA.9i . we 
obtain a set of self-consistent equations for GJJg.(a;). In the 
wide band limit, Taai'^) = Tq, we get 



a;±iO+ -uo' 



■dw'. 



(A. 16) 



Similarly, we have I]2^°(a;) — ±.\TPp°' [lu). Collecting all 
these expressions we finally obtain 



GS.(c.) = 



w-£d,-ReI]is(cjff)+ir(l + /(cjff))-2irPj:(LJs) 

(A.17) 
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